Nonlinear microscopic relaxation of uniform magnetization precession 



Vladimir L. Safonov and H. Neal Bertram 

Center for Magnetic Recording Research, 
University of California - San Diego, 
9500 Gilman Drive, La Jolla, CA 92093-0401 
(February 1, 2008) 



Dynamic relaxation for nonlinear magnetization excitation 
is analyzed. For direct processes, such as magnon-electron 
scattering and two-magnon scattering, the relaxation rate is 
determined from the linear case simply by utilizing the magne- 
tization oscillation frequency for nonlinear excitation. For an 
indirect process, such as slow-relaxing impurities, the analysis 
gives an additional relaxation term proportional to the exci- 
tation level. In all cases the effective magnetization damping 
is increased compared to Landau-Lifshitz-Gilbert damping. 



I. INTRODUCTION 

The study of damped magnetization dynamics of fine 
ferromagnetic particles and thin films is important for the 
development of nanomagnetic devices and high-density 
magnetic recording. A conventional theoretical tool to 
study magnetization relaxation is based on the phe- 
nomenological Landau-Lifshitz equation jjj or its mod- 
ification with Gilbert damping These equations 
conserve the absolute value of magnetization (|M| = 
const) in a single domain region due to a strong ex- 
change interaction approximation. They are relatively 
simple and therefore have been utilized for various cal- 
culations and micromagnetic simulations. However both 
Landau-Lifshitz and Gilbert equations were introduced 
(a) for small magnetization motions and (b) for the case 
of high magnetic symmetry (axial symmetry) with an 
isotropic damping fitting parameter a ("damping con- 
stant"). Nevertheless both "a" and "b" conditions are 
usually violated. 

Recently a theoretical approach || , Q , H , Q , has 
been developed to correct the limitations of the Landau- 
Lifshitz-Gilbert (LLG) theory. The main idea was to 
represent the magnetization dynamics as the motion of 
a damped nonlinear oscillator with the random force of 
thermal fluctuations. The oscillator model is a conve- 
nient tool to establish a "bridge" between the microscopic 
physics, where the rotational oscillator variables a* and 
a naturally describe spin excitations (as creation and an- 
nihilation operators) and the macroscopic magnetization 
dynamics. In particular, it has been rigorously shown 
by including specific coupling of a magnetic system to 
a variety of loss mechanisms, that for small oscillations 
near equilibrium the macroscopic damping term reflects 



the anisotropy of the system || , ^ . Only for complete 
uniaxial symmetry do the LLG results apply. 

Nonlincarity in the relaxation process appears with 
the increase of magnetization deviation from equilibrium. 
Depending on the situation, the nonlinear magnetization 
damping either increases or decreases the total relax- 
ation rate. For example, in our simulations of magne- 
tization reversal (large magnetization motion) byspin- 
wave dynamics in a quasi-single-domain grain |10| we 
have demonstrated that an effective damping parameter 
a can be about hundred times greater than that in the 
linear regime. For a single-domain grain a phenomeno- 
logical extension of the oscillator approach to nonlinear 
excitations has been made . In this case for large am- 
plitudes the effective damping is much greater than that 
of LLG. 

The aim of this paper is to give a microscopic analysis 
of nonlinear magnetization damping. We consider both 
direct and indirect loss mechanisms and show explicitly 
how nonlinearity is included in the relaxation rate. The 
approach is to first transform the magnetization dynam- 
ics without damping to rotational oscillator coordinates. 
This transformation yields equations where the oscilla- 
tor frequency depends on the degree of nonlinearity. The 
analysis parallels the approach for low-level excitations. 



II. THE NONLINEAR OSCILLATOR MODEL 

Here we review the magnetization dynamics in terms 
of a nonlinear oscillator. For simplicity we shall con- 
sider a single-domain grain in the case when the external 
magnetic field Ho is parallel to the uniaxial anisotropy 
axis (z). To describe the magnetization dynamics M(i), 
we introduce the classical spin S = M.V/hj, where V 
is the grain volume, U is Planck's constant and 7 is the 
gyromagnetic ratio. Using the Holstein-Primakoff trans- 
formation: 



S + = a V2S - N, S z 



S~ = a* V2S - N, S d 



S — a* a, 
= S x ±iS y , 



(1) 



where N — a*a — \a\ 2 (N < 25), we can represent the 
magnetization dynamics in terms of the oscillator vari- 
ables: 



da/dt = —iuj(N)a. 



(2) 



The oscillator (g) is characterized by the magnetization 
oscillation frequency for nonlinear excitation: 



1 



u(N) = wo[l - N/S(l + Ho/H K )}, 



(3) 



which reflects the change of effective magnetic field 
with increasing magnetization deviation from equilib- 
rium. Hk is the anisotropy field, ujq — j(Hq + Hk) 
is the ferromagnetic resonance frequency. 

Previously it was assumed that the corresponding 
stochastic differential equation is of the following form 



(d/dt + T)(N)) a = -~iu{N)a + f(t), 



(4) 



where f(t) describes a random uncorrelated force and 
ri(N) is the nonlinear damping. The formulas (|I]) and 
(0) are valid within one potential well up to the top 
of the energy barrier < N < N top , where uj(N) = 
(N top — 5(1 + Hq/Hk)Y Utilizing the fluctuation- 



dissipation theorem for ( 
given by |J: 

t](N) = ati(N), 



the nonlinear damping is 



a 



»/(0)/wo. 



(5) 



III. NONLINEAR RELAXATION MECHANISMS 

In this section we show how nonlinear relaxation arises 
from microscopic mechanisms. The linear relaxation rate 
77(0) is well understood and determined by solving the 
coupled magnetic - thermal bath equations. The main 
mathematical idea is that the principal complex am- 
plitude dynamics in the small damping approximation 
(r](N) <C w(iV)) is given (from either (||) and (^)) by 



a(t) oc exp[— iui(N)i\. 



(6) 



For all cases where the relaxation is of a direct form of 
coupling (see, ||), the microscopic analysis will yield the 
same relaxation expression for a given mechanism except 
that luq is replaced by oj(N). 



A. Magnon-electron scattering 

Let us consider the magnon-electron scattering process 
in a ferromag netic metal §, Q. A magnon with wave 
vector k = and energy Tiuiq and a conduction electron 
with wave vector k ^ and energy Tilo^ are transformed 
into a conduction electron with wave vector k' and en- 
ergy hujek' — hLu e k + fiijQ. This confluence process occurs 
in the presence of defects to violate momentum conser- 
vation. The linear relaxation rate for this process can be 
written as: 



T]m — e 



k,k' 



/kk' 



N 



(n k - n k i)6(uj ekl - cj ek - uj ), 



(7) 



where Cdef is the defect concentration, /kk' /N is the scat- 
tering amplitude and 7Tk is the Fermi occupation num- 
ber. Neglecting ojq dependence in the scattering am- 
plitude and taking into account that ujq -C lu^ and 
nk — TTk' — —uodrik/dujek, from (^) we obtain: 



Vm-e = Cdef^o 



duj e ir 



E 



kk' 



N 



^ek) 



(8) 



with a linear dependence 7/„ 
u>(N), we have: 



oc ujq. Replacing ujq 



f] m -e(N)/uj(N) = r/ m - e /u = a m _ e (9) 
in agreement with the phenomenological relation (^) . 

B. Two-magnon scattering 

The linear regime of two-magnon scattering on defects 
has been considered by many authors (e.g., |l^] , jl3| , , 
|]l5|). The linear relaxation rate for this process is given 
by: 



V2r, 



= 77 E l G k| 2 <5(w k - Wq), 



(10) 



where Gk describes the scattering amplitude and Wk is 
the spin- wave frequency with k ^0. The abovemcntioned 
simplification gives: 



m m (N)=TvJ2\Gk\ 2 S[oj k (N)-ui{N)} 



(11) 



Here u>k(N) is the spin- wave frequency taking account 
of nonlinear excitation. This example does not give, in 
general, a linear relation between r]2 m (N) and ui(N). 



C. Slow relaxation 

The energy loss in this case occurs via an intermediate 
damped dynamic system, the 'slow relaxing' impurities 
(e.g., ||, p3| l, Jl6[). The magnetization motion modu- 
lates the impurity splittings and varies the thermal equi- 
librium populations of the energy levels. Let us consider 
two-level impurities with energy: 



(12) 



where Jloj is the splitting frequency, rij is the upper lever 
population and j is the impurity index. The impurity 
level modulation is defined as: 



(13) 



2 



The nonlinear dynamic equation for the complex am- 
plitude is: 



da/dt = —iuj(N)a — id(7ii mp j /Ti)/da* 
= -iuj{N)a - i^^drijit). 



(14) 



The kinetics of the impurity population is defined by the 
following equation: 

drij/dt = — Tuj(Qj)[rij — tit (CI j)]. (15) 

Here riij(fij) is the impurity relaxation rate and 
riT(Clj) = [exp(hflj i 'ksT) + is the equilibrium pop- 
ulation at frequency Clj (t) = fioj + dClj (t) . We can we 
solve Eq.dTl) and obtain: 



5rij(t) 



eMhdfr) ~ hAV) rnj[n 3 -(*i)] (ie) 



x{n T [Clj(ti)] - n T (O ,j)}dii, 
where Srij(t) — nj(t) — tit(CIoj) and (j>\\,j(t) = 
J T^Q^dt'. 

— oo 

In the case of small modulation ft\8Vlj(t)\/k-B,T <C 1 we 
can write the following expansion: tit^Iq j + SClj(t)] = 
nr(fio,j) + [dn,T(Cloj) /dClo,j]Sflj(t) + ... Neglecting for 
simplicity the Qj(t) dependence of Fuj, we can rewrite 
Eq.(|l6[) in terms of a(t) and a*(t) and substitute these 
terms into d||). Thus, we can obtain the nonlinear re- 
laxation rate in the form: 



where 



2 ( dn T {Clo,j)\ ^o^\\,j 



3 



is a linear relaxation rate and 



'Is 



(17) 



(18) 



(19) 



scattering, the nonlinear relaxation rate was shown to 
be directly proportional to the nonlinear frequency as in 
the initial phenomenological approach (0). For the case 
of two-magnon scattering the relaxation rate, in general, 
is not proportional to the nonlinear frequency. A spe- 
cific case of indirect interaction, 'slow-relaxing' impuri- 
ties, was analyzed in detail. The resulting relaxation rate 
was the sum of the linear term plus a nonlinear term pro- 
portional to the level of excitation. An absence of linear 
relation between the nonlinear damping 77(A) and mag- 
netization oscillation frequency for nonlinear excitation 
u(N) indicates a colored thermal noise in the system. 

In Ref. [§ it was shown that the nonlinear oscillator 
damping is greater than the conventional damping in the 
Landau-Lifshitz equation. This result is a general con- 
clusion for all the relaxation mechanisms considered here. 
Corresponding to (|j) and (^) the magnetization dynamic 
equation can be written as 

dM/dt = -7M x H eff + (a/M s )M x dM/dt, (20) 



where, in general, 



N M,-M z 



u)(N)(l- N/2S)' S 



M, 



1 — cos ( 



(21) 



Here 6 is the deviation angle. In the vicinity of equilib- 
rium (A = 0, M z = M s , 6 = 0) a = a = 7/(0) /a; , where 
a is the LLG parameter. Away from equilibrium, the 
damping parameter increases 5 with increasing nonlin- 
earity, becoming infinite as N — > 2S (M z — ► — M s , or 
8 — > tt). An application of this approach to magnetiza- 
tion reversal is discussed in Ref. p7|. 
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is the coefficient of a nonlinear damping. We can see that, 
even though the relaxation rate depends on the level of 
excitation, in general, r] Br (N)/uj(N) ^ r} sr /uio- 



IV. DISCUSSION 

We have considered nonlinear relaxation for both di- 
rect and indirect coupling to a thermal bath. For di- 
rect coupling it was shown that the linear relaxation rate 
can be simply converted to the nonlinear rate by replac- 
ing the linear frequency ojq by the nonlinear frequency 
ui(N). In one case of direct coupling, magnon-electron 
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